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ABSTRACT 



Supersonic combustion in the inviscid stream has been 
shown to reduce the base drag on a projectile. Fast reacting 
fuel burned externally in the inviscid, supersonic flow gen- 
erates compression v/aves, which, when properly originated, 
impinge on the recirculation bubble aft of the projectile 
base and increase the base pressure. A free- jet wind tunnel 
was constructed to study external burning where heat addition 
was simulated by nozzle contour. Sufficient increase in base 
pressure was attained to not only reduce base drag but also 
provide thrust. Thus, external burning provides a method 
for increasing the performance of projectiles. 
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I. 



INTRODUCTION 



For an axisyiranetric object flying at supersonic speed as 
much as 50 percent of the overall drag can be attributed to 
the negative gage pressure acting upon the base. This rela- 
tionship between base drag and the remaining drag on the body, 
which consists of wave, skin friction, and boat-tail drags, 
is observed primarily in the low supersonic flight regime. 

For Mach numbers near 2 the pressure on the base can be re- 
duced to as little as 50 percent that of the ambient air 
while the nose pressure can be kept relatively low by using 
slender, streamlined shapes. With increasing Mach number, 
however, a limit to the base drag is reached since the base 
pressure cannot be less than a negative one atmosphere, 
whereas the pressure on the nose v;ill continue to increase 
with Mach niomber. 

The reduction of drag is desirable in order to increase 
the range and terminal velocity of a projectile. Pressure 
or wave drag and skin friction drag have been satisfactorily 
reduced by designing streamlined shapes. Base drag, however, 
has not been reduced significantly although the use of boat- 
tails and base bleed have contributed slightly. The most 
promising method of base drag alleviation is that due to 
external burning. Heat addition in the inviscid flow in 
proximity of the base will create compression waves that, 
if optimally originated, will impinge on the base flow 
region. The interaction of these compression waves will 



11 



deflect the free streamline toward the symmetric axis, and 
the static pressure in the recirculation zone will increase 
thus reducing the base drag. 
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II. THEORY 



The base flow problem is visualized most easily with 
the inviscid flow model in Fig. 1 which ignores the effects 
of viscosity. 

A. INVISCID FLOW 

A semi-infinite, two-dimensional model can be formu- 
lated which easily satisfies all the necessary boundary 
conditions . 

1 . Tv70-Dimensional Model 

If a value of (= base pressure) is selected, 
then the initial turning through the Prandtl Meyer expansion 
is uniquely determined. The pressure and, therefore, the 




FIG. I 

TWO-DIMENSIONAL INVISCID FLOW 
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velocity and Mach number are constant along the streamline 
BC. The Mach number can be calculated from the Prandtl 
Meyer equations. For uniform, two-dimensional flow over a 
convex corner, the pressure depends only on the angle of 
inclination of the free streamline (angle 0 in Fig. 1) , and 
thus BC is a straight line. The shock wave originating at 
C is weakened eventually by the interaction of the expansion 
waves. Near the plane of syimnetry, the shock must return the 
flow to the free stream direction. As the flow continues 
downstream of the trailing shock, the static pressure ap- 
proaches the free stream pressure. Unfortunately, for a 
given M^ , an infinite number of possible solutions exist. 

As P„ is decreased, the Prandtl Meyer angle is increased, 
and point C moves closer to the base. A limit is reached 
eventually where for a given local Mach number M^, the 
trailing shock cannot return the flow to its original 
direction. The maximum deflection shock wave for a given 
local Mach number determines a limiting flow from the infi- 
nite number of possible solutions. This limiting case 
yields the minimum possible P_ and thus the maximxam base 
drag in inviscid flow. 

2. Axisymmetric Model 

Inviscid theory can be extended to axisymmetric 
shapes with some complications. The expansion characteris- 
tics issuing from the base corner are not straight lines as 
in the Prandtl Meyer flow except near the corner. Because 
the flow conditions upstream of the trailing shock are 
dependent not only on the streamline inclination at the 
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base corner but also on the history of the flow upstream of 
the Mach lines passing through the corner, the free stream- 
line of constant pressure is not straight, as shown in Fig. 
2. The method of characteristics can be utilized in order 
to determine possible flow patterns. 




The inviscid flow pattern of an axisymmetric body cannot be 
constructed to the axis of symmetry and still satisfy the 
boundary conditions. This is a consequence of the curvature 
of the free streamline and the fact that for a given Mach 
number, a single trailing shock is limited in the amount of 



15 



turning it can provide. The angle of inclination of the 
constant pressure free streamline continues to increase 
monotonically as the axis of symmetry is approached. The 
flow pattern, therefore, cannot be continued beyond a point 
where the limiting angle (a in Fig. 2) exceeds the angle 
that the single trailing shock wave can return the flow 
parallel to the axis. In real flows, the wake prevents the 
free streamline from reaching the axis, and thus this limita- 
tion yields some similarity to observed flows. Flow fields 
containing a free streamline which does not reach the axis 
of symmetry can be considered as those that would exist in 
an inviscid flow about an axisymmetric body which has an 
infinitely long cylindrical sting attached at the base as 
shown in Fig. 3. 
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The fraction d/2h equals the ratio of rod diameter to 
body diameter. As d/2h approaches one, the tv70-dimensional 
result for is obtained. With a rod present, inviscid 
theory predicts a lower base pressure and consequently a 
higher base drag for two-dimensional profiles than for 
axisymmetric shapes. 

This is in qualitative agreement with the experiment. 
Chapman [Ref. 1] found that the base drag of axisymmetric 
bodies in an inviscid flov; was so small that it could not 
account for observed drag values determined from experi- 
ment in real fluid flov;. This strongly suggested that 
viscous effects were essential in determining base flow. 

B. VISCOUS FLOW 

The inclusion of viscosity effects results in a more 
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causes a slow circulatory motion in the dead air region, and 
the viscous mixing process causes the boundary layer to 
thicken as it approaches the trailing shock wave. The base 
pressure arrives at a steady state equilibrium value because 
of the scavenging effect of the outside flow on the recircula- 
tion region. As some of the dead air mass is removed, the 
angle of turning at the base corner is increased, and con- 
sequently the base pressure is decreased. The larger turning 
angle increases the flov; velocity outside the boundary layer 
which increases the scavenging action, thereby lowering the 
base pressure. The process is repeated until point C in 
Fig. 4 reaches an equilibrium position. As point C moves 
toward the base, however, the pressure ratio of the trailing 
shock wave increases, making it more difficult for the 
scavenged air and the low velocity air in the boundary layer 
to overcome the pressure rise. An equilibrium position is 
reached v;hich balances these two opposing effects. Viscosity, 
therefore, determines a unique solution from the infinite 
number possible for a given Mach number and Reynolds number 
but, unfortunately, complicates the theory excessively. 

C . CHAPMAN ' S CORRELATION 

Viscous theory is too involved to be a practical 
method for estimating base pressure. Instead many correla- 
tions have been attempted to aid in the understanding or to 
create a practical means of predicting base drag. According 
to Chapman's model [Ref. 1] base pressure depends on the 
Mach number and pressure at the base corner , the boundary 
layer thickness <f, the body diameter 2h, and the body angle 
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^ at the base. For finite bodies, as shown in Fig. 5, 
profile shape directly affects the Mach number and pressure 
at the base corner since the existence of a bow shock wave, 
for example, would alter the Mach number and pressure prior 
to the corner. 




M, and P, are needed instead of M and P . Chapman 

11 oo oo 

offered a simple means of calculating the corrected approach 
Mach number M^ and pressure P^. He suggested extending the 
after body one diameter (2h) and taking the average Mach 
number and pressure along this extension as the corrected 
Mj^ and P^. From dimensional analysis. Chapman found that 
c B ) where C^g is the base pressure coeffi- 
cient and is equal to (Pg-P^j^) / ( /J vi/^2 )• Viscous effects 
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were included in a Reynolds nuinber dependency of the boundary 
layer thickness S. Reynolds number. Re = /^Vo 5 >L//ta, , is based 

upon the finite body length L and the free stream conditions. 

is the free stream air density, Vo© is the velocity, and 
Aft? is the coefficient of viscosity. Assuming a fixed Mq,, 
boundary layer theory predicts boundary layer thickness S is 
proportional to l/'AfRe for laminar flow or S is proportional 
to l/([Re)^^ turbulent flow. With yS= 0, the final result 

from dimensional analysis is ^ laminar flow 

and C turbulent flow. 

Chapman used data from axisymmetric flow to make a corre- 
lation of Cpg with L/KRe''^ for laminar flow, as shown in 
Fig. 6, and L/KRe"^^) for turbulent flow. 

The effect of both Re and body shape is apparent, and the 
laminar correlation is reasonably close to the data. For 
turbulent flow. Chapman's correlation again is good; however. 

Re has a much smaller effect on base pressure. This is con- 
sistent v/ith theoretical considerations since is much 

less Re dependent than is 
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FIG. 6 

MEASURED AND CORRELATED BASE PRESSURE 
DATA FOR M^=l.53 LAMINAR BOUNDARY 
LAYER FLOW 
REF I 
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D. EFFECTS OF SPIN 

The effect of spin, or a body's revolution about its 
axis of symmetry, on base pressure also has been shown to be 
considerable. Unclassified material on this subject is 
lacking, and apparently no purely theoretical approach has 
ever been made successfully. It has been concluded instead 
of analysis, that in order to determine the effects of spin, 
experimental measurements on specific object shapes need to 
be made . 

C. Wieselsberger [Ref. 9] showed in 1927 that the drag 
coefficient of bodies at subsonic speeds could be increased 
by spin. Eased on drag data of a projectile and cylinder at 
low subsonic speeds v/ith circumferential velocities up to 
three times the fom-zard velocity, Wieselsberger concluded that 
the increase of total drag resulted from increases in skin 
friction and base drag due to spin induced alterations of the 
boundary layer velocity profile. R. Lehnert [Ref. 9] showed 
that spin influenced the aerodynamic characteristics of pro- 
jectiles at supersonic speeds, but comprehensive information 
concerning the effect of spin was lacking. H. H. Kurzweg 
[Ref. 7] showed a close relationship between boundary layer 
and base pressure, as did Chapman, which prompted a study by 
R. Lehnert and S. Hastings [Ref. 9] to determine the base 
pressure behavior as a functional of spin rate. Experiments 
were conducted at the Naval Ordnance Laboratory in a super- 
sonic aeroballistics tunnel. Cone cylinders of varying 
lengths were spun at rates varying from 0 to 43,000 rpm, and 
base pressure was measured through an orifice located on the 
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sting immediately adjacent to the model base. Spark -schlieren 
photographs were taken in order to study the effect of spin 
rate upon the wake . 

In Fig. 7 base pressure ratio v/as plotted versus model 
spin rate with the dimensionless model length as a parameter. 




The upper curve, v;hich corresponds to the shortest model, 
shows a considerable decrease in base pressure ratio with in- 
creasing spin rate. The next lower curve also shows a decrease 
in Pg/P^ with increasing spin rate. The three lower curves, 
corresponding to the three longest models, show relatively 
little change in the base pressure ratio with changing spin 
rates. As was determined by both Kurzweg and Chapman, for a 
given xMach number and model configuration, the shape of the 
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curve for base pressure versus Reynolds number is determined 
by the state of the boundary layer. 

Figure 8 shows base pressure ratio plotted against 
Reynolds number for Lehnert and Hastings' data. 




The Reynolds number was referred to free stream conditions 
and the wetted model length. The first curve shows the zero 
spin data. The large decrease in with increasing Re up to 
3.6 million is accompanied by the existence of a laminar 
boundary layer along the entire model length. Further in- 
creases in Re resulted in a decrease in P , and the transition 

B 

in the wake mixing zone appeared to move upstream. Base 
pressure reached its minimum value when transition occurred 
at the model base and increased only slightly at higher Re 



24 



when transition occurred upstream of the base. The curves 
for various spin rates show a considerable decrease in base 
pressure with increased spin when laminar flow existed at 
the base. This decrease in base pressure due to spin v/as 
accompanied by an upstream travel of transition in the mixing 
zone as was the case for zero spin. In the transition region 
the base pressure is unaffected by spin when the boundary 
layer is turbulent at the base. Increasing spin resulted 
in a slight base pressure increase. The minima of the base 
pressure versus Re curves move toward lov/er Re V'/ith increas- 
ing spin rate. It can be concluded that spin destabilizes 
the boundary layer and thereby produces premature transition. 
Transition at zero spin was observed to occur at nearly one 
diameter upstream of the base while at 40,000 rpm, the transi 
tion moved to almost two diameters upstream. 

Kurzweg [Ref. 8] and Chapman [Ref. 1] observed that base 
pressure and wake length were correlated. Decreasing base 
pressure corresponds to decreasing wake length as shov'/n in 
Fig. 9. This vras confirm.ed by Lehnert and Hastings [Ref. 9] 
and suggested a method of measuring base pressure on free 
flight projectiles. Note the similarities of Fig. 9 with 
Fig. 7. 
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WAKE length vs. model SPIN RATE 

REF 9 



E, HEAT ADDITION 

Base drag reduction by external combustion v/as studied 
by VI. C. Strahle [Ref. 13]. Strahle observed that the injec- 
tion of mass or heat into the viscous recirculation bubble 
resulted in low efficiency since the base pressure rise was 
limited to the free stream atmospheric limit. Combustion in 
the inviscid supersonic stream, however, imposed no such 
limit and, therefore, showed the most promise for reducing 
the base drag or for providing a propulsive mechanism. 

1 . Control Volume Analysis 

A control volume analysis [Ref. 6] as shown in 
Fig. 10 gives some intuitive feeling for heat addition in 
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the supersonic stream. Applying the momentum theorem for 
the control volume in the x direction, one has 

dA - j[V«*-c/A P- A (0 

where dA is the projected area perpendicular to the free 
stream direction x. For planar flov; and no heat addition, 
pressure P along streamline cd is determined by the local 
flow angle in accordance with the Prandtl Meyer function. 




Figure 11 shows that the static pressure has a lower value 
in the fast turning case (short recirculation region) as com- 
pared to the slower turning case (long recirculation region) . 

The will be less for a fast turning stream- 

c 

line; and, consequently , P., also will be lower. Since by the 
continuity equation fjuiA = f%iA , only changes in momentum 
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thickness £** and displacement thickness can change the 
momentum flux through de and be, it is the gage pressure 
integrated over cd and de that largely determines base 




Compression waves emitted from the heat addition zone as 
shown in Fig. 12 increase the static pressure along the edge 
of the shear layer, cd in Fig. 11; and therefore, rcp-R^idA 
is increased thereby increasing Pg. The pressure decreases 
along cd as the flow turns toward the axis. By adding heat 
in the inviscid region, the shear layer is turned inward 
with nd or less loss of static pressure along cd. Pg is in- 
creased which results in a base drag reduction. 
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2 . Strahle Two-Dimensional Model 

Strahle [Ref. 13] restricted his heat addition 
model to two-dimensional, turbulent, supersonic flow, with 
the heat addition zone aligned with the inviscid stream 
beyond the base corner turn as shov;n in Fig. 13. 

The heat addition zone is rectangular in shape and is 
assumed far enough removed from the viscous region to pre- 
clude any heated streamlines being entrained in the viscous 
flow region, which would result in a performance loss. 
Strahle chose the mixing theory of Crocco and Lees [Ref. 2] 
to describe the base flow because of the method's success in 
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predicting the base pressure. Strahle found that the 
inclusion of external combustion did not alter the base 
flow theory of Crocco and Lees but instead only changed the 
initial conditions for Crocco - Lees governing differential 
equation (Eq. 2.30 in [Ref. 2]). 




With heat addition, the angle of inclination O in Fig. 

14 is no longer a function of alone but is also a function 

of both W and x. 
e 

The symbol W^ is the ratio of the x component of velocity 
to the speed of sound a^. A saddle point singularity exists 
through which the trajectory must pass in order that the far 
wake limit can be reached. This forms the uniqueness 
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condition; since without it, an infinite number of solutions 
would exist. This saddle point is called a critical point 
and physically represents the transition point from subsonic 
to supersonic flow in the viscous region. A unique solution 
is forced by the critical point regardless of the external 
combustion occurring in the inviscid stream. The position of 
the critical point, however, and the influence on the allowable 
initial condition will be altered since & is influenced by 
combustion . 




Strahle analyzed the combustion region by the linear theory 
of heat addition applied to a supersonic stream. The X and Y 
coordinates defined in Fig. 14 are obtained by rotating the 
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axis through an angle 0(o)and are aligned v;ith the rectangular 
heat source. In the uniform field v/ith Mach nxrmber equal to 
M_ , the perturbation velocity potential equation is 

''Pq 0VY “-f - C2) 

and the pressure coefficient is = “2.^>c/Ua. Upon inte- 
grating Eq. (2) where f is non zero only for positive X and 
Y, along the Y = 0 line, the perturbation velocity potential 

1 = 0y()(,o) = - l//3lfJfCK C3) 

With 0(^0") unknown and thus the orientation of , which is the 
heat release per unit volume per unit time, also unknown, the 
Bernoulli equation along a streamline states 

<lP/fa=-ciu/u = = -<I(M')/2M'‘Q+U-0kV2] (4) 

Solving forJ© and integrating yields 

e-6t = (5) 

where 0 is the Prandtl Meyer function. 0 ^ is evaluated at 
0 ( 6 ) on the Y=0 line regardless of the history of the viscous 
region edge. For large heat addition, this is a source of 
error since the initial direction of the viscous region is 
different from its alignment at the conclusion of the in- 
fluence of heat addition. Integrating Eq. (3) and substi- 
tuting into Eq. (5) gives 



e-6>b- 






where n. = ( 8 -0 W 




gCx") = O 


if 


X< X, ^ X > 


= X-X;-^gY; 


if 


X, ^ X < Xj 


= /^8 ^ 


,-f 


X < X 3 


“ CV; - X 4- X; Wj 


if 


^3 ^ ^ 
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with 



X, = X, ^ 

Xj. = Xi + + 

Xj^Xi+W-f^jY; 

^4 “ ^ i t W ■*■ |5 j5 ( V; + 



3 . Strahle One-Dimensional Model 

Because the linear theory contains contradictions 
of the assumptions at even modest levels of heat addition, 
Strahle suggested a one-dimensional approach for comparison. 
Figure 15 shows the basic scheme under a small angle assiimp- 
tion that =©a.-^i and dAc/d>« 

v/here 7v. is the area of the heat addition zone and A is the 
h ^ 

area of the adiabatic free stream. 







FIG. 15 

strahle's one-dimensional keat addition 
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Assuming the flow on the upper side of (Fig. 15) 

contains waves of only one family, then ©j, 4 ’ . 

Using the one-dimensional flow relations that CdAe/Ae')cl = 
0-Me)?rK&Pa/Pc')cj>i , Strahle obtains 

0-©t,=^O(Mb)~ C>a^V(:Ae/PG)[0'{^e")AMc]JPe/cly>-<^Aj^/dK . (?) 

dPe/dx is always maintained small by the combustion-affected 
region. Thus, for large heat release, dominates the 

pressure gradient term. If is identified with dAj^/dx 

in Eq. (6) , then an analogy exists between the one-dimensional 
and linear theory approaches. On the average, 
b(AV/v)/W where is the constant pressure volume expan- 

sion of the combustible mixture, b is the heat addition zone 
capture height, and W is the length of the combustion zone. 
With assumed in the strongest region v;here X, ^ X ^ Xj, 

in Eq. (6) , the equivalence is made that 

AV/V = ~ 0 AK^//< 0 % and a = A V/vW/9g C8) 

where A is the heat release per unit mass of fuel, ^ 

is the air to fuel ratio, and . Functionally, 

tv/v = f(A, Pe.Te ) with the most important dependences 
being on M. and Tg, which is the base temperature. Strahle 
hoped that by combining linear theory with a one-dimensional 
approach the features of both methods could be included v;hile 
removing most of the objections of each method individually. 

4 . Strahle *s Performance Parameters 

Strahle computed the force on the base to be 
Fe OCPa-PeJK where P^ is the pressure on the base (where 
X = 0) and h is the base half height. The specific impulse 
is defined as where is the fuel flow rate. 
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For small values of force level on the base and letting /t->co 

I.P' C2AhfC?-i)/M8^„tW](P3/P,)^JP/dn) (9) 

where P = Pg/Pj^ = a heating parameter. The fact that 

/{-►«> creates an absolute linear heat addition limit, which 
though impractical , does suggest the strong role that the 
critical point plays in the nonlinear regime. 

Figure 16 shov;s the effect the heat addition zone location 
has on specific impulse. As is moved from the base toward 




STRAHLES SMALL DISTURBANCE SPECIFIC IMPULSE df^d^ 

VS. COMBUSTION ZONE LOCATION X,- 



the critical point, which is located at approximately x = 7.7 
in the linear limit, the specific impulse (iP/An) reaches a 
maximum when the combustion zone is roughly centered between 
the base and the critical point. A steep drop in effective- 
ness occurs as the zone nears the critical point and extends 
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beyond it. This is expected because any influence upon the 
reattachment process dov;nstream of the critical point where 
the flow is supersonic cannot be transmitted upstream. The 
linear limit predicts that external combustion will result 
in a base pressure rise, the specific impulse is expected to 
be high, and the combustion is wasted if the length of the heat 
addition zone W is too long. 

Strahle recomputed the example in Fig. 16 for finite 
combustion strength, and the resulting axial pressure dis- 
tributions are shown in Fig. 17. 




FIG. 17 

STRAHLE's FINITEAXIAL PRESSURE DISTRIBUTIONS 



Four of the curves are for = 2.0 which corresponds to 
the maximum I location in Fig. 16. The top curve is for 

= 0.2. Apparently, for this case, it is best to have X^ 
small. The reason for this lies in the nonlinear behavior 
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of the critical point. In all cases, the critical point be- 
comes entrained in the expansion portion of the combustion- 
affected zone. The combustion length W, therefore, fixes the 
interaction length, and Strahle suspected that the most effi- 
cient placement of the combustion region v7ould be as close to 
the base as possible. He also reasoned that the base pressure 
increase resulted from the entrapment of the subsonic recircula 
tion bubble by compression v/aves emitted by the combustion zone 
The singularity can occur only v;hen the inviscid flow is 
directed inward or for 0^0 . The viscous zone thickness, 
hov;ever, must remain finite. Since the primary effect of 
the combustion region is the downv/ard turning of the flow, the 
last available position for the critical point appears to be 
within the combustion-affected portion of the flow. The 
flow passing the base corner cannot turn as far without com- 
bustion, v;hich in turn indicates a rise in base pressure. 

5 . Axisymmetric Heat Addition 

Strahle 's two-dimensional heat addition requires 
modification for axisymmetric application. A method by Fuhs 
[Ref. 6] suggests that the change in base drag resulting from 
external burning is v;here A„ is the area of the 

base for both planar and axisymjnetric bodies as shovm in 
Fig. 18. 

Equal amounts of heat are released in both the planar 
and axisymmetric cases. Heat addition in the planar case is 
represented by the lines above and belov; the base area and 
in the axisymmetric case by a concentric circle of radius r^. 
Since in both cases, the length of heat addition zone is 
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, the base area influenced by heat release is in the 

ratio 



axisymmetric base area .. rrr rr r r 

planar base area ' / 2- • ^ /2 



CIO) 




FIG. 18 

COMPARISON! OF BASE! ARFA FOR PLANAR AND 

AXI SYMMETRIC GEOMETRY 



The axisymmetric geometry apparently has less base area 
influenced by heat addition. Flow through an annular heat 
addition region v/as investigated. It v/as found that the 
compression waves merged sooner v;ith axial symmetry, and that 
in contrast with the planar case, there is an optimal r 2 /rj^ 
ratio. 
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III. APPARATUS 



In order to study the effects of heat addition on the base 
pressure of an axisyminetric body, a free jet wind tunnel was 
constructed at the Naval Postgraduate School's blow down 
laboratory. The test model simulated a 5 -inch projectile at 
23,000 feet and a Mach number of 2.0, which are typical condi- 
tions for a maximum range trajectory. Similarity parameters 
were matched by using a Mach 2.0 free jet coaxial nozzle 
shown in Fig. 19 which gave a Reynolds number based on radius 
of 1.51 X 10®. 




FIG. 19 

COAXIAL CYLINOER/NOZZLE 
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A turbulent boundary layer v/as assvuned, and the displace- 
ment thickness^* was computed to be 0.041 inches at the 
projectile base. This provides a similarity parameter of 
displacement thickness to radius ($*/h of 0.0164. 

The equipment consists basically of a 300 psig air supply, 
a globe valve, a pressure reducer valve, a 90° elbow, a center 
cylinder support/flow straightening section, a nozzle and an 
overhead exhaust hood. A 4 -inch pipe feeds the air from the 
supply tank to the on/off globe valve. The pressure reducer 
valve is followed by a 90° elbow which directs the air flov^? 
vertically and provides a feedback pressure for the regula- 
tion of the pressure reducer valve. A 2 -inch spacer is 
located between the elbow and the center cylinder support 
section and provides internal access to the centerline cylinder. 
The centerline cylinder is supported inside the center cylinder 
support section by flow straighteners located at the top and 
bottom of the 3 -foot section of 4 -inch pipe. A section of 6- 
inch pipe houses the nozzle on top of the cylinder support 
section and provides armor protection and a support for the 
pressure probe as shown in Fig. 20. Figure 21 shov;s the 
exhaust hood above the apparatus, the mercury manometer board 
used to obtain pressure data, the mercury vapor light source 
and one mirror used in the schlieren system for flov7 
visualization . 

Figure 22 shows the cylinder extending from the cylinder 
support section with the nozzle and armor section removed, 
along with the pressure gage connected to the air supply 
reservoir . 
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FIG. 21 

PHOTOGRAPH OF TEST APPARATUS 
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ri6.22 

TEST CYLINDER 



43 



The cylinder base is formed by a cap which contains six 
pressure taps. Four taps are located across the base radius, 
and two taps are located just upstream of the base corner. 

The pressure from each tap is carried dovmward through 
the cylinder and the cylinder support in plastic tubes v;hich 
exit the apparatus through the 2 -inch spacer located on top of 
the elbow. The tubes are connected to a mercury manometer from 
which the pressure data was read. Total and static pressure 
also can be measured by the external pressure probe, which 
will be the primary base pressure source when cylinder spin 
is included. The results of this study were for the static 
case with all similarity parameters matched except for spin. 
However, because spin appears to have a considerable impor- 
tance on base pressure, provisions v/ere included to allow 
spinning the projectile. When a 5-inch projectile is fired, 
the barrel rifling spins the projectile to 250 revolutions 
per second. This establishes a spin similarity parameter of 
r<^/u = 0.1595 where r is the radius of the projectile, cO is 
the angular velocity and u. is the free stream velocity at 
Mach 2.0 at 23,000 feet altitude. The similarity parameter 
will be matched in the free jet experiment by spinning the 
cylinder at 48,000 rpm. The spinning cylinder design con- 
sists of a 6.5-inch stainless steel cylinder which slides 
over a support cylinder containing high pressure nitrogen 
orifices to provide gas bearing support and align the 
spinning cylinder radially. A small internal impulse turbine 
driven by high pressure nitrogen will spin the cylinder. 
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Heat addition v;as simulated by nozzle contour. As de- 
scribed in Ref. 6, the external burning process can be 
simulated by nozzle wall contours v/hich generate character- 
istics identical to the supersonic combustion of an external 
burning assisted projectile. A gas generator within the 
projectile injects fuel-rich, hot, combustible mixtures into 
the supersonic inviscid flow v;here ignition occurs due to 
high temperature and the chemical nature of the fuel. The 
region of heat addition was assumed axisymnietric about the 
cylindrical body of the projectile with the radius in 
Fig. 18 equal to the body diameter. For stoichiometric fuel 
mixtures, calculations shov;ed that the combustion would result 
in a 4.5° deflection of the streamline bounding the heat 
addition zone, and it was this 4.5° deflection inward of the 
streamline at T 2 from the cylinder centerline which governed 
the wall contour. Half stoichiometric nozzles also are 
planned for future analysis. An iterative procedure using 
the method of characteristics was used to design the nozzles 
to give the proper -4.5° flow angle. Using V7 to denote the 
length of the combustion zone (see Fig. 14) and h, the cylinder 
radius, several nozzles of varying V7/h v/ere fabricated with an 
epoxy cast poured over a wax plug to the correct contour. 

These nozzles provide several data points shown in Table I 
which show the effect external heating has on base pressure. 
During testing, spacers were used to vary the axial location 
of the heat addition zone with respect to the base corner. 
Figure 23 is a photograph of the nozzles, spacers, and tunnel 
armor section. 
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N. W/h 

VW 


0.72 


1.84 


4,34 


-2.23 


0.703 


0.6 


1.145 


-0.95 


0.6 85 


0.858 


1.253 


-0.5 


0.759 


0.96 


1.20 


0.0 


0.8S3 


1.055 


1.14 


0.5 


0.993 


1.074 


1.005 




TABLE I 
TEST MATRIX 
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FIG. 23 

PHOTOGRAPH OF NOZZLES, SPACERS, AND TUNNEL ARMOR SECTION 



The symbol 1 represents the distance from the corner 
where the first fully developed compression wave impinges 
the cylinder v/all or the wake downstream of the base, and h 
again is the cylinder radius as shown at the top of Figs. 

40, 41, and 42. Five 1/h values were observed for each 
no 2 zle as shown in Table I. Negative 1/h denotes compression 
beginning upstream of the base corner, positive 1/h signifies 
compression beginning downstream of the corner, and 0 denotes 
compression starting at the corner. Table I is the test 
matrix obtained to date showing the results of the completed 
tests . 
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IV. TEST PROCEDURE 



After the proper nozzle and spacer combination v;as 
installed on the apparatus and the air reservoir was charged 
to 300 psig, a data point for the test matrix was taken. The 
procedure involved opening the globe valve and setting the 
proper output pressure from the pressure reducer valve, A 
pressure tap in the elbow v/as connected to a separate 
pressure gage which was used to set the conditions for each 
run. This tap also provided feedback pressure for the 
pressure reducer valve. With the proper output pressure set, 
a schlieren concave mirror system was used to photograph 
each test run as shovm in Figs. 25 - 39 . Polaroid positive/ 
negative film was used in order to facilitate print 
reproduction. Base pressure information was read from a 
mercury manometer and recorded for data reduction. At the 
conclusion of the data readings, the pressure reducer valve 
was turned off, and the air reservoir was recharged for the 
next run. 
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V. RESULTS AND CONCLUSIONS 



The pressure ratios (P„/P_ ) as shown in TalDle I indicate 

jj oo 

the performance of the observed simulated external burning. 
Figure 24 is a plot of the results. 




The data points taken for the two uniform flow nozzles 

compare favorably with M = 2 results compiled by Przirembel 

and Page [Ref. 10] . These points represent the base pressure 

ratio of the cylinder with no heat addition; and because 

P_/P is less than 1.0, base drag results. All other data 
B oo 

points, except one, show an increase in base pressure ratio 
with the applied combustion, indicating a reduction in base 
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drag. A13. the data for W/h = 4,34, and the last tv/o points 
for W/h = 1.84 in fact, indicate thrust developed at the base 
since Pj^/P^ for these cases is greater than 1.0. The data 
indicate, as expected, that a longer heat addition zone 
(Larger W) will result in a higher Pj^/Pqo ij^^i^^ase, since a 
longer interaction length of compression v/aves acting upon 
the recirculation bubble has more effect than a shorter zone. 

It also appears that maximum benefit from the heat addition 
occurs v/hen the entire recirculation bubble is embedded in 
the combustion-created compression v;aves. For V7/h = 4.34, 
a decline in performance is noted as 1/h exceeds -1.0, which 
indicates that for a long heat addition zone, optimal axial 
placement is necessary so as to not v;aste compression v;aves 
downstream of the sonic point. Because supersonic flow exists 
dovmstream of the sonic point, any compression effect occurring 
beyond that point will not be transmitted upstream into the 
subsonic circulation bubble, which v;ould be necessary for a 
base pressure increase. For the shorter heat addition zones 
(W/h = 1.84 and 0.72), the base pressure ratio increases as 
the zone is moved rearward. It is assumed that no compression 
was wasted beyond the sonic point up to 1/h = 0.5. Although 
not tested, a decline in performance is expected v/ith further 
rearward movement of the heat addition region v/ith respect to 
the cylinder base corner. Figures 25-39 show a schlieren- 
derived photograph of each data point obtained with a 
pictor’ial description of the heat distribution actually 
simulated. The double cross-hatched regions depict anomalous 
compression, and the single cross-hatched areas show the heat 
addition zones without multiple reflection. 
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FIG. 25 

PHOTOGRAPH AMD SCHEMATIC FOR w/h- 4.^4, -2.' 
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F\G. 26 

PHOTOGRAPH AND SCHEMATIC FOR W/h= 4.34^ i/k= -0.9S 
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ri G. 27 

PHOTOGRAPH AND SCHEMATIC PC R W/K= 4.34^ J?/K= -0.50 



54 




FIG. 28 

PHOTOGRAPH AND SCHEMATIC FOR W/K » 4 . 34 , = 0.0 
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FIG. 29 

PHOTOGRAPH AND SCHEMATIC FOR W/K= 4.34^ VK- O. SO 
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fig. 30 

photograph AND SCHEhAATIC oF W/f^= 1.84^ J?A = -2. 23 



57 




n G. 31 

PHOTOGRAPH AND SCHEMATIC FOR W/K^ 1.84, M = -0.95 
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riG. 32 

PHOTOGRAPH AND SCMETMATlC FOR W/K = 1. 04, je/M-0.5o 
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FIG. 33 

PHOTOGRAPH AMO SCHEMATIC FOR W/^ = I. 84-^ J?/h = 0.0 



60 






n G. 34 

PHOTOGRAPM AND SCHEMATIC FOR W/K= 1.84, Jf/h = 0.50 



61 



I 




-rrpi^ 

I 





riG. 35 

PHOToeRAPM AND SCHEMATIC FOR W/K = 0.72, V^> = "2.23 
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FIG. 36 

PHOTOGRAPH AND SCHEMATIC FOR W/h J?/h = -0.9S 






-^= 0.759 

Ro 




FIG. 37 

PHOTOGRAPH AND SCHEMATIC FOR W/h =Q 72^ JP/h = -O.SQ 
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FIG. 38 

PHOTOGRAPH AND SCHEIMATIC FOR W/h =0.72, JP/K= O. O 





fig. 39 

PHOTOGRAPH AND SCHEMATIC FOR W/l^ =0.72^ O. SO 
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A Study of the pictures shows that a uniformly distrib- 
uted heat source v^7as not obtained as desired. Instead, the 
starting compression waves, which originated immediately 
downstream of the nozzle throat, reflected from the cylinder 
v;all to the nozzle corner and then into the base flow region 
as shownat the top of Figs. 40, 41, and 42. For the W/h = 
4.34 nozzle, the simulated heating distribution is shown in 
the middle of Fig. 40. 

Compression, originating from the simulated heat addition 
zone, starts at point b and increases to uniform strength at 
point e where it levels at a constant valvie until just prior 
to point h. The reflected compression v;aves from the nozzle 
corner superimposed upon the uniform heat cause a spike which 
is approximately twice the magnitude of the intended uniform 
value. This spike is seen in the W/h = 4.34 photographs as 
a shock wave, formed by a coalescing of the reflected 
compression waves, as seen in Figs. 25-29. 

Figures 41 and 42 show the expected compression wave re- 
flection pattern and heat distribution for W/h = 1.84 and 
W/h = 0.72 respectively. Basically, the same phenomena of 
compression wave reflection occurred. The shorter heat addi- 
tion zones, however, ended prior to the reflection, and, 
therefore, the reflected waves created an almost equal 
strength combustion region downstream of the intended heat 
addition zone rather than a spike as seen for the W/h = 4.34 
case . 
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The photographs shov/ that for the vVh = 1.84 case, the 
reflected compression waves never coalesce into a shock but 
are spread over a slightly longer axial distance. For W/h = 
0.72, a shock did form, but its strength is notably less than 
the shock occurring for W/h = 4.34. 

Although the reflected compression waves existed for each 
data point, it is assumed that the base pressure was affected 
by the reflection for only the points circled in Fig. 24. 

The circled points were determined by the location where the 
reflected wave impinged the recirculation bubble. As pre- 
viously described, base pressure can be influenced by 
compression waves only if the waves interact with the base 
flow prior to the critical point. Any interaction dovaistream 
of the critical point will not influence pressure at the base. 
The bottom of Figs. 40, 41, and 42 shows the estimated 
critical point location for the three W/h nozzles. 

A vertical line extended through the intersection of a 
given 1/h value and the critical point curve will have on 
its left the portion of the heat addition region affecting 
the base pressure. 

In conclusion, it is evident that external burning is 
capable of not only reducing base drag but also producing 
thrust under the proper conditions. The optimal benefit is 
derived when tlie entire recirculation bubble is embedded in 
combustion created compression v/aves. Waves beginning up- 
stream of the base corner have limited value for raising base 
pressure, and waves interacting dovmstream of the critical 
point are wasted. 
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